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RADIATION FROM SLOTS ON CYLINDRICAL BODIES 
USING GEOMETRICAL THEORY OF DIFFRACTION 
AND CREEPING WAVE THEORY* 

By Constantine A. Balanis 
Langley Research Center 

SUMMARY 

A hybrid solution employing wedge diffraction and creeping wave theories is used 
to compute the radiation patterns of axial and circumferential slots in the principal planes 
(equatorial and elevation) on conducting cylinders of finite and infinite lengths. The slots 
are excited by parallel -plate waveguides operating in the TEM and TE^q modes. 

For the equatorial -plane pattern, the total field in the lit region is obtained by the 
superposition of two fields, that is, the wedge -diffracted and the creeping-wave fields. 

The wedge-diffracted field is obtained by approximating the parallel -plate — cylinder 
geometry with two wedges, each formed by a wall of the waveguide and a tangent plane to 
the cylinder surface at the edge point. The creeping -wave contribution is obtained by the 
use of diffraction and propagation coefficients of waves traveling around conducting 
curved surfaces. The total field in the shadow region is obtained solely from the 
creeping-wave contribution. For the elevation-plane pattern, wedge diffraction tech- 
niques for the entire pattern are employed. 

The method is checked computationally by comparison with the modal solutions for 
axial and circumferential slots on right circular cylinders and experimentally for ellip- 
tical cylinders since modal solutions are not readily available. Experimental models 
are also used for the verification of the elevation plane pattern computations since 
boundary -value solutions are not available for finite length cylinders. 

Computed results using this technique compare favorably with those obtained from 
existing modal expansion boundary -value solutions and experimental results. The main 
advantages of the present technique are that it can be applied to geometries where modal 
solutions are not possible, in numerical ranges where the convergence properties of 
modal expansions are relatively poor, in parametric design problems since the contribu- 
tion from each field is separated, and in the analysis of antennas with finite physical 
sizes. 

*The material presented herein was offered as a dissertation in partial fulfillment 
of the requirements for the Degree of Doctor of Philosophy, The Ohio State University, 
Columbus, Ohio, June 1969. 



INTRODUCTION 


Since the geometrical theory of diffraction has been introduced, it has been 
employed successfully in the solution of various types of diffraction problems. Previous 
work in applying the wedge diffraction and creeping wave theories, each one individually, 
to scattering and radiation antenna problems is discussed. No one has attempted a 
hybrid solution with both wedge diffraction and creeping wave theories used for the same 
scatterer. The work presented in this report is an attempt to formulate a hybrid solu- 
tion employing both wedge diffraction and creeping wave theories on the same antenna. 
The antenna is a slotted, axially and circumferentially, conducting cylinder, with the 
slots being excited by parallel -plate waveguides operating in the TEM and TE^q modes. 

The radiation properties of slot antennas on conducting cylinders have been studied 
extensively. (See refs. 1 to 5.) Such analyses have been restricted to geometries where 
modal expansion field functions are readily available and whose physical length is infinite. 
Moreover, the solutions obtained exhibit poor convergent properties for objects more 
than a few wavelengths in extent. However, many practical problems involve slot 
antennas mounted on finite-size conducting bodies whose geometrical shape does not con- 
form to a coordinate system where the wave equation is separable. Wedge diffraction 
(refs. 6 to 9) and creeping wave theories are used in the analysis presented herein to 
compute the radiation patterns of axial and circumferential slots on conducting cylinders 
of finite and infinite lengths. The method employed is an extension of the geometrical 
theory of diffraction (refs. 10 and 11) and creeping wave theory (refs. 12 to 16). 

The geometrical theory of diffraction and creeping wave theory are extensions of 
geometrical optics which account for diffraction. They introduce diffracted rays in addi- 
tion to the usual rays of geometrical optics. These diffracted rays are produced by 
incident rays which hit edges, corners, vertices, or curved surfaces. Some of these dif- 
fracted rays enter the shadow regions and account for the field intensity there. The dis- 
tribution of the incident energy among those diffracted rays is described by a diffraction 
coefficient. The diffraction of the incident rays is a local phenomenon, so the diffraction 
coefficient depends on the geometry of the edge or curved surface in the vicinity of the 
point of diffraction and on the polarization of the incident field. 

A field is associated with each diffracted ray and the total field at a point is the 
sum of the fields of all rays at that point. The phase of the field on a ray is assumed to 
be proportional to the optical length of the ray from some reference point. Appropriate 
phase jumps must be added as a ray passes through a caustic (congruence of rays at a 
point or line). The amplitude is assumed to vary in accordance with the principle of 
conservation of energy in a narrow tube of rays. The initial value of the field on a dif- 
fracted ray is determined from the incident field with the aid of an appropriate diffraction 
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coefficient. These diffraction coefficients are determined from certain canonical 
problems. 

The canonical problem, which yields the diffraction coefficient for an edge, is the 
diffraction of a linearly polarized plane wave by a two-dimensional, infinite wedge. The 
diffraction coefficient is found from Sommerfeld's asymptotic, high frequency solution to 
this problem. (See ref. 17.) Pauli (ref. 18) has obtained a practical formulation for the 
diffraction coefficient and Hutchins (ref. 19) has generalized it. He has obtained a series 
solution which is valid for the exterior and interior regions of a wedge. In addition, his 
solution is accurate for wedges of large included angles whereas the Pauli solution is not 
accurate. Oberhettinger (ref. 20) has obtained a similar series in which the leading 
term is identical to the Fresnel integral form for half -plane diffraction. The work of 
Hutchins tends to bridge the solutions of Pauli and Oberhettinger. 

The primary task in applying the creeping wave theory is the determination of the 
diffraction, attenuation, and ray-path factors for a general body (and it is not yet always 
practical). Thus, it is necessary to evaluate these factors for canonical targets whose 
exact solutions are available, such as the cylinder and sphere. (See refs. 12 to 16 and 
21 to 23.) In order to obtain a more general solution for these factors it is convenient 
to utilize experimental data to validate an empirical model for more general targets such 
as the prolate spheroid (refs. 24 to 26). 

The attenuation and diffraction coefficients for a cylinder have been obtained by 
comparing the asymptotic series of the scattering boundary -value solution with that of 
the diffraction problem formulation (refs. 13 to 16). The ray paths along the surface of 
the cylinder are geodesics (great circles). The scattering by a sphere is, however, a 
most important problem. The classical solution of Mie using separation of variables 
and series techniques can be considered as the starting point. Senior and Goodrich 
(ref. 21) have obtained an asymptotic form to the Mie series solution for the sphere 
through the application of the Watson transformation. Hong (ref. 22) has derived attenu- 
ation and diffraction coefficients with higher -order correction terms. The formidable 
task of implementing this solution numerically has yet to be accomplished. 

An approach to the general formulation of the creeping -wave paths on bodies of 
revolution has been developed by Kinber (ref. 27). Kinber has shown that the wave equa- 
tion may be expanded in a set of ray coordinates in which the solution can be written in a 
form whose magnitude is dependent upon the cross section of the ray tube and whose 
phase is dependent only upon the propagation constant and path length traversed. 

The empirical approach utilizes a simplified ray -path geometry to determine 
approximate attenuation and diffraction coefficients for a sphere. (See refs. 24 to 26.) 
This approach is of interest in that it lends itself for extension to more general targets, 
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and is a simplified creeping wave analysis. It uses a single nonconvergent (and non- 
divergent) ray path together with approximate diffraction and attenuation coefficients to 
construct an approximate solution for the scattered field due to creeping waves propa- 
gating in each direction along the ray path. The path chosen is the path traversed by the 
TT major M ray (that is, the path corresponding to the E-plane of the sphere). This analysis 
suggests that an approximate picture of scattering by a sphere can be constructed by 
neglecting the creeping waves which have a radial magnetic field (that is, the "minor" 
creeping waves) and by considering only the creeping waves which have a radial electric 
field (that is, the "major" creeping waves) except for the H-plane scattered fields. 

The prolate spheroid represents a body for which an exact closed -form solution is 
not available except on the axis of rotation. In order to find the ray-path geometries of 
the prolate spheriod for an arbitrary angle of incidence, the geodesic corresponding to 
the point of attachment and the tangent direction at that point must be calculated. The 
determination of the tangent direction of the creeping wave at the shadow boundary is not 
an easy task in general. It is therefore suggested that the simplified ray-path geometry 
(refs. 24 to 26) for the sphere can be applied to the prolate spheriod. In this case, the 
creeping -wave paths are elliptical; thus, the total attenuation must be expressed as an 
integral which is dependent upon the radius of curvature along the path. Also the radius 
of curvature at the points of attachment and reradiation must be computed in order to 
determine the diffraction coefficient. The specific solution for the prolate spheriod has 
been presented by Ryan (refs. 24 and 26) and Peters and Ryan (ref. 25). 

Many other antenna problems have been treated using wedge diffraction techniques. 
Russo, Rudduck, and Peters (ref. 28) applied the geometrical theory of diffraction to cal- 
culate the total antenna pattern of a horn in the E -plane, including the backlobe region. 
Because the theoretical and experimental patterns are in excellent agreement, the method 
for treating diffraction by edges is valid. Obha (ref. 29) used the geometrical method of 
diffraction to calculate the radiation pattern and gain of a finite width corner reflector 
antenna. This method also yields good results for the computation of backscattering 
from an antenna having conducting plates finite in extent. 

Extensive work of applying the wedge diffraction method to waveguides has been 
carried out by Rudduck and coworkers (refs. 6 to 9). The principal tool employed is 
diffraction by a conducting wedge; the resulting electromagnetic field may be treated as 
a superposition of the geometrical optics field and the diffracted field which behaves as 
a cylindrical wave radiating from the edge of the wedges. These techniques may be 
applied to any two-dimensional antenna or scattering body which may be constructed 
from a set of wedges; the radiation or scattering pattern may be obtained for any excita- 
tion which can be expressed in terms of plane or cylindrical waves by superposition of 
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the individual wedge diffractions. Some elementary examples include parallel -plate 
waveguides, walls of finite thickness, and polygonal cylinders. 

A basic feature of this technique is that it does not require an assumption of the 
value of the total field or current on some surface, as do conventional methods of aper- 
ture integration which employ the approximation of physical optics. This feature thus 
allows more accurate treatment of problems than ordinarily obtained by approximate 
methods; it also provides knowledge of the fields in terms of the incident field, which 
results in the ability to analytically determine admittance and gain of antennas (ref. 6), 
radiation patterns (refs. 7 and 8), mutual coupling (refs. 30 and 31), and reflection coef- 
ficient (refs. 32 to 34). 

The geometrical ray techniques of this method provide conceptual simplicity with 
which solutions may be formulated. This, together with superposition of wedges, allows 
structural aspects to be taken into account; for example, antennas mounted with and with- 
out ground planes, structures with thin or thick walls, and arbitrary waveguide trunca- 
tions can be treated. In addition, the fields in all space can be determined; thus, back- 
lobe regions can be treated. 

The essential feature of this approach is that it is applicable to all types of com- 
plicated problems, some of which are not tractable. The approach is to resolve a 
complicated problem into simpler ones, each of which will have a relatively simple solu- 
tion. The basic building blocks will be rigorous solutions to such canonical diffraction 
problems as the wedge, cylinder, and sphere, which can be put together to solve more 
complicated problems. 

The work outlined in this hybrid solution is for general conducting cylinders of 
arbitrary convex cross section. However, computations are restricted to circular cyl- 
inders for which boundary -value solutions exist for comparison and elliptical cylinders 
for which models and experimental data are available for comparison. The radiation 
patterns of interest for each case are the principal -plane patterns, equatorial plane and 
elevation plane in the usual spherical coordinate system. 

The equatorial -plane pattern is obtained by the superposition of diffracted fields 
from a set of wedges and creeping -wave fields propagating around the surface of the 
cylinder. The set of wedges approximate the parallel -plate — cylinder geometry in the 
vicinity of the discontinuity formed by the waveguide and the cylinder. Each wedge is 
formed by a wall of the parallel -plate waveguide and the tangent plane to the cylinder 
surface at the edge point. 

For the equatorial -plane pattern, two distinct regions can be identified: the lit 
and shadow regions. The shadow region encompasses the space occupied by the set of 
wedges and cylinder surface and the lit region the remaining space. The total field in 
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the lit region is obtained by the superposition of two fields - that is, the wedge-diffracted 
and creeping -wave fields. The total field in the shadow region is obtained solely from 
the creeping -wave contribution. For the elevation -plane pattern, wedge diffraction tech- 
niques for the entire pattern are employed. 

The most interesting part of this hybrid solution is the coupling mechanism between 
the wedge -diffracted and creeping -wave fields. Creeping waves are launched by the tip 
of the wedges in a direction tangent to the cylinder surface at the edge point. Since the 
field must be continuous along the line separating the lit and shadow regions, the wedge- 
diffracted fields along the surface of each wedge will be used as the initial value of the 
creeping waves to preserve the continuity of the fields. The creeping waves continually 
radiate energy in a tangential direction as they propagate around the cylinder surface 
until they reach the opposite wall of the slot and again illuminate the wedges; this is a 
second-order field contribution and it is neglected. The loss of energy due to reradia- 
tion is accounted for by the use of an attenuation factor. 

The method is checked computationally by comparison with the modal solutions for 
axial and circumferential slots on right circular cylinders and experimentally for ellip- 
tical cylinders since modal solutions are not readily available. Experimental models 
are used for the verification of the elevation -plane -pattern computations since boundary- 
value solutions are not available for finite-length cylinders. 

Computed results from this technique compare favorably with those obtained from 
existing modal expansion boundary -value solutions and experimental results. The main 
advantages of the present technique are that it can be applied (1) to geometries where 
modal solutions are not possible, (2) in numerical ranges where the convergence prop- 
erties of modal expansions are relatively poor, (3) in parametric design problems since 
the contribution from each field is separated, and (4) in the analysis of antennas with 
finite physical sizes. 


SYMBOLS 

No attempt has been made to define all the symbols used in this report in a single 
list; most symbols used in this analysis are defined by diagrams for the various geom- 
etries. The following symbols are used as main symbols with appropriate subscripts 
throughout the analysis: 

D diffraction coefficient 

E electromagnetic field 


6 



creeping -wave field 


wedge -diffracted field 
geometrical optics field 
total field 
diffracted ray 
reflected ray 
attenuation constant 

phase constant 
wedge angle 

Numerical subscripts indicate particular wedge. 

Superscripts within parentheses indicate order of diffraction. 

WEDGE DIFFRACTION 

Wedge diffraction techniques are applied for the analysis of the radiation pattern 
of the parallel -plate waveguide shown in figure 1. This geometry has been treated by 
using the wedge diffraction concepts by Rudduck and his coworkers (refs. 6 to 9). Their 
solution is incorporated in the slotted-cylinder radiation-pattern analysis and is outlined 
in this section. A single -diffracted wave emanates from each wedge that is illuminated 
by an incident plane wave. The single -diffracted waves produced by one wedge illuminate 
the other with double -diffracted fields being produced. This process continues to higher 
orders of diffraction. 


E, 


Et 


Erp 

r d 

^R.FL 

a 


k 

WA 


Single Diffraction 

The principal method employed in the analysis of a parallel-plate waveguide is dif- 
fraction by a conducting wedge. The diffraction of a plane wave by a wedge was solved 
by Sommerfeld (ref. 17). Pauli (ref. 18) obtained a practical formulation of the solution 
for a finite angle conducting wedge which was later improved by Hutchins (ref. 19). (See 
the appendix.) The total electromagnetic field from the wedge may be treated as the 
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Figure L- Geometry of a parallel-plate waveguide aperture. 





Figure 2.- Geometry for plane wave wedge diffraction with geometrical optics 

region. 
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superposition of the geometrical optics and the diffracted fields which behave as cylin- 
drical waves radiating from the edge of the wedges. 

The diffraction of a plane wave by a wedge is shown in figure 2. The solution to 
the plane-wave diffraction problem may be expressed in terms of scalar functions that 
represent the normal and parallel polarization components of the electromagnetic field 
to the plane of study in figure 2. The total field is defined as 

E = Eq + Ejj (1) 

where Eq is the geometrical optics field and Eq is the diffracted field. The dif- 
fracted field is given by 

E D = V B ( r » ^~^o> n ) ± V B( r ’ ^o> n ) < 2 ) 

where the parameters r, \(/, \p Q , n are shown in figure 2 and Vg(r, ty±\p 0 , n) is 
the diffraction function as defined in the appendix. The plus sign applies for the polari- 
zation of the electric field normal to the edge 

^ =0 (3) 

9n 'wedge 

where n is the unit normal to the edge of the wedge. The minus sign applies for the 
polarization parallel to the edge 


E wedge ® 


(4) 


The geometrical optics field is defined in three regions. (See fig. 2.) For plane-wave 
incidence, the geometrical optics field is 


Geometrical _ 
optics field 

Incident Reflected 

geometrical field ± geometrical field 

(Regions) 


Eq 

e jkrcos(^-i// 0 ) 

(IT - l// Q < ^ < 77 + l^ 0 ) 

(5) 

Eq 

e jkrcos(i/'-i/' 0 ) ± e jkrcos(i//4-i^ 0 ) 

(0 < 1// < 7T - l^o) 

(6) 


e g = 0 

The time dependence ej w t is used throughout this analysis. 


(jr + ) (V 
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The diffracted wave Ej) may be represented as a cylindrical wave radiating from 
the edge. (See the appendix.) In fact, at large distances from the edge and in regions 
removed from shadow boundaries, Ep has the radial dependence e“j* cr / \fr. Because 
of this cylindrical nature, subsequent diffractions of a diffracted wave may be treated as 
the diffraction of a cylindrical wave by a wedge. 

Diffraction by Pair of Wedges 

The process of diffraction by a pair of wedges (parallel -plate waveguide) is now 
outlined. Two polarizations which must be distinguished are those corresponding to the 
TEM and TE^q waveguide modes shown in figure 3. In the TEM mode, the incident plane 
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(b) TEio mode. 

Figure 3.- TEM and TEiq modes in a parallel-plate waveguide aperture. 



wave is parallel to the axis of the guide having a polarization perpendicular to the edge of 
each wedge. The TE-^q waveguide mode may be represented by two TEM waves which 
reflect obliquely back and forth between the waveguide walls. (See ref. 35.) Thus, in the 
TEjq mode the incident plane waves have a polarization parallel to the edge of the wedges. 

TEM mode.- In the TEM mode an incident plane wave propagates parallel to the 
axis of the guide walls as shown in figure 3(a). The far -zone singly diffracted fields 
from wedges 1 and 2 are given by 


E Dl( r l^l) = 


Ef e 


A slni 

ni n. 




27rkr 1 


v + 


cos 


ni 


cos 


( 8 ) 


E D2 ( r 2 ,< ^ > 2) ~ 


E 2 e 


i-t 


\/ 27rkr 2 


COS 


n. 


- COS 


7 T - (pc 


(9) 


The superscripts denote the order of diffraction. If the phase of the incident waves with 
respect to edge CD is expressed as 


E x = E = 1 


( 10 ) 


and 


Eg = Ee 


-jkwcotc^ -jkwcot<£ 

& e ® 


then equations (8) and (9) reduce to 

E I)l( r l>*l) 


- j ( kr l + j) 


R 


( 1 ) 


^2vrkr 1 D1 


(*l) 


E D2 ( r 2’ ^ 2 ) 






R 


27rkr r 


( 1 ) 

D2 


(*2) 


( 11 ) 


( 12 ) 


(13) 
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where R^ (0jJ and (0 2 ), the sin S 1 y diffracted rays from wedges 1 and 2, are 

given by 


,< 1 ). 


R 


( 1 ) 

D1 


(♦l)" 


-J-sin-1- 
n l n l 


IT 77 + ^1 

COS - — cos 


n l 


ni 


( 14 ) 


^2 (*2) = 


-jkw cot 0 g j 
e 5 — sin 71 


n 2 n 2 


7T - 

COS — - cos 

n 2 n 2 


(15) 


Applying the far -field approximations 


r l “ r o + 2 cos (0o + 0g) 
r 2 * r o - | cos (0 O + 0g) 


(16) 


for phase terms and 


r r r 2 K r o 


(17) 


for amplitude terms with 


01 * 02 * 00 


(18) 


gives 


E Dl( r o-*o) - 


l{k[r o +j cos (*o+0g|+jj 

/2vkr^ 




(19) 
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and 


E D2 ( r o’ *o) _ 


{k[r 0 + |cos( V ^ g jj + f 


^27rkr c 


<20) 


The singly diffracted rays may again be diffracted to produce doubly diffracted rays and 
so on to higher orders of diffraction as shown in figure 4. The singly diffracted ray 

illuminating edge (2) gives rise to the doubly diffracted ray R^; R^ causes R^ 

in a similar manner. Also, some of the singly diffracted rays from edge (I) are reflected 


( 2 ) 

r di (( V 



Figure 4.- Doubly diffracted and reflected rays of a general parallel-plate waveguide. 
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from wedge 2 and appear to radiate from the image of edge (I) ; these reflected rays are 
given as 


R RFL (* 0 ) ~ R D1 ( _< M (fg < $0 < |) ( 21 ) 

The doubly diffracted ray from edge © is expressed as 

R D2 (*o) = «DlG|yB( h . V* V n 2) + V B( h ' »- n 2)] ( 22 > 

where 

R D1G = R D1 ("^g) ( 23 ^ 

and from edge (I) as 

4 2 1 (*o) = 4 1 2 G |yB( h ' *o + *g. ” 1 ) + V B( h > ”l)] (24) 

where 

R D2G = R D2( ir-< ^ > g) ( 25 ) 

An additional doubly diffracted ray by edge (I) from the image source is given by 

R d!r K) = Hgipf B ( 2w > t *o> n l) + V B ( 2w > T + * 0 ’ n l)] < 26 ) 

where 

e dip = r di(-|) < 27 ) 

The total diffracted rays from edges (D and ©, by using single-double diffractions, are 
expressed as 

R DlK) = “Sw + 4 2 l(*o) + ^r(*o) (28) 
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( 29 ) 


and the reflected ray from wedge 2 as 

R RFl(^o) = R Dl("^o) 

The total diffracted field from the aperture may be expressed as the superposition 
of the total diffracted rays from edges (D and (2) plus the total rays from the image 
source to yield 


-jfk[>o + f cos ('V 4 t')) + j] r 

= 5 The; F“ ( *°> 

jkhcos(0 o +^ K ) -j2kw sin <p3 

+ r D 2 W<* v e * R DiHo) e J (31) 


Each term in equation (31) contributes to the radiation pattern only in certain regions as 
follows: 


Ray 

Region 

R Dl(^o) 

-<Pg < <p Q < 7T - WA1 

R D2 (Oo) 

IT + WA2 < <j> 0 < TT - (f> 

r D1 Ho) 

+ d> < d) < — 
^ g v o 2 


There are subsequent diffractions which result in third- and higher-order diffrac- 
tions from edges ® and (2) . The total higher-order diffractions (that is, second and 
higher order) can be put in a closed form. (See ref. 6.) The total illumination of edge © 
from edge Q can be expressed as 

R iG“ R iHe) <32) 

where Ri(0 o ) is the total diffracted ray from edge (J). Consequently, the total higher- 
order diffractions from edge (2) are given by 
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R D2 CM = R lo|yB( h ' ^o^g’ n 2) + V B( h ’ 7r “ W n 2) 


(33) 


and the total ray by 

B 2(M' R $(*o) + B 82(*o) < 34 > 

The total higher -order illumination of edge (I) is given by 

R 2G =R 2 (’ r -'M (35) 

and 

r ip = r i(-|) < 36 > 

Thus, the total higher -order diffractions from edge 0 are given by 

R D1 K) = ^G^B^’ n l) + V B( h ’ 27r+ ^o"^g’ n l)J 

+ r ip|v b (2w, |+ 0 O , n^ + V B (2w, |^+ 4 > 0 , n x jj (37) 

and the total diffractions by 

R i(^o) = R mVo) + R Di('i , °) < 38 > 

The total diffracted field from the aperture can be expressed as the superposition 
of the total diffracted rays from edges 0 and 0 plus the total reflected ray as given by 
equation (31) where (^o)’ R D2 (0o)> and R Dl( _< ^o) are replaced by R i (0 O )> 
r 2 (0 o )> and R i( _< ^o)> respectively. 

The total diffracted waves from edges 0 and 0 are given in terms of the unknown 
illuminating rays Rjq, R 2G 3 and R 1P* These rays can be determined by the solution 
of three simultaneous linear equations formed by expressing each unknown ray in terms 
of equation (34) or equation (38) and are given as follows: 

R 1G = *S}g + R 2G V 2G (-*g) + R 1P V 1P (-*g) < 39 > 
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(40) 


R 1P ” ^lP + R 2G V 2G (' 2) + R 1P V 1P (‘ 2 ) 

R 2G = R D2G + R lG V lG( 7r-< ^g) 

where the quantities Vjq, V^p, and V^q are the unit-wave diffractions used in equa- 
tions (33) and (37). 

TEjq mode.- The diffraction at the aperture of the parallel -plate waveguide for the 

TE 10 mode may be treated in a similar manner as the TEM mode. The TE^q mode may 
be represented by two plane waves reflecting obliquely back and forth between the wave- 
guide walls at an angle 


A 0 = sin -1 (42) 

2w 

as shown in figure 3(b). The wave has a polarization parallel to the edges of the wedges 
which form the waveguide walls. For this mode two cases must be distinguished as 
follows: 


Case I: A 0 > 4>g 

Case II: A 0 S (pg 


For case I, edge (2) is not illuminated by the incident plane wave and no singly diffracted 
ray emanates. For case II, both edges are illuminated by the incident wave; hence, 
singly diffracted rays emanate from both edges. 

Since the polarization is parallel to the edges of the guide for the TE-^q mode, the 
negative sign in the diffraction formula is chosen. The singly diffracted ray from edge (J) 
is obtained from equation (2) as 


R mW “iq-slni.. 


cos 


- cos 


+ 4>o - A o 


ni 


cos 


ni 


- cos 


7T+ 


+ 


A 


(43) 


ni 
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For case n (A 0 S (p\ f the singly diffracted ray from edge © is given by 


41(V) = -sr sin £ 


* -jkw(sinA 0 + cot 0gCOsA o ) 

1 sin — e ' ' 


/ 1 

1 cos — - 
\ n 2 

IT - 



n 2 


TT - 


COS — - COS 
n 2 


Po * A o 


( 44 ) 


where the exponential factor and minus sign represent the phase of the incident plane 
wave at edge (J) . The reflected rays are given by 


R RFL if o) - " R Dl(-0o) < ^o < |) ( 45 ) 


where the preceding minus sign results from the reflection. 

Multiple diffractions occur in the same manner as for the TEM mode but with the 
minus sign chosen in equations (22), (24), (26), (33), and (37). Thus the total higher-order 
diffracted waves from edges (J) and (2) for the TE^q mode are given by 


R D1 M ~ R 2G [ V B( h ’ ^o +< ^g’ n l) “ V B (h, 2 7r+ (pQ-ipg, n^jj 
- R 1 p[ V B ( 2w ’ f + n l) - V B (2w, |4 (pQf n x | 

and 


(46) 


R D2 (■#>«) 



*-4> 0 -<l> g , n 2 ) 



w-^o +< ^g» 



(47) 


The minus sign preceding R^p results in the same manner as for the reflected rays of 
equation (45). The total wave from each edge is obtained by using the TE^q mode rays 
in the same equations valid for the TEM mode, that is, equations (34) and (38) to (41). 
The unknown illuminating rays are determined in the same manner as for the TEM mode 

by using the formulations for Rp], r d 2 > R d 1» and R D2 S iven for the TE^q mode. 


CREEPING-WAVE DIFFRACTION CONCEPT 


When a wave is incident upon an opaque object which is large compared to the wave- 
length, a shadow is formed. However, some radiation penetrates into the shadow region 
due to diffracted rays as shown in figure 5. (See refs. 12 to 16.) These rays are pro- 
duced by incident rays which are tangent to the surface of the body. Each tangent ray 


18 




Figure 5.- General concept of creeping-wave diffraction. 

splits at the point of tangency with one part continuing along the path of the incident ray 
and the other traveling along a geodesic on the surface of the body. At each following 
point, it splits again with one part traveling along the geodesic and the other reradiating 
along a tangent to the geodesic. From a single incident ray, infinitely many diffracted 
rays are produced, one of which is reradiated at each point of the geodesic. These waves 
traveling around the opaque body have been designated as creeping waves introduced first 
by Franz and Deppermann (ref. 12) for the interpretation of scalar diffraction by circular 
cylinders and spheres. 

The scattered field caused by the creeping wave mechanism for a plane wave inci- 
dent on a cylinder is given by (ref. 13) 


-jk(t+s) v V^ ds 


E c (s ,<f>) = Ej(Q) 5 ^ D mh (Q) D mh (P) e 
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which for a right circular cylinder reduces to 


E c (s,0) = Ei(Q) 


e -jk(t+s) 

w~ 



( 49 ) 


since the radius of curvature is constant. The diffraction coefficients and attenuation 
constant are given by 


a]L 

12 


D mh = D mh (Q) D m h( p ) = Diffraction coefficients » 

= Attenuation constant ~ (kp) (4m + N)J e"* ® 


-(kp) 1/3 d m 

(50) 

2/3 jf 
e 6 

(51) 


where 

p radius of curvature of body 

k phase constant 

t path length along body 

s distance from point of detachment to observation point 
P point of detachment of creeping wave 

Q point of attachment of creeping wave 

For a hard surface (E-field normal to the surface) d 0 = 1.083, d^ = 0.555, and N = 1; 
for a soft surface (E-field parallel to the surface) d Q = 0.645, dj = 0.490, and N = 3. 
The diffraction coefficients and attenuation constants are obtained by comparing the 
asymptotic series expansion for large values of kp of the canonical boundary -value 
problem with equation (48). (See refs. 13 and 14.) 

Rays diffracted by the edges of a slotted cylinder (shown in fig. 6) which are tan- 
gent to the surface of the cylinder behave in a similar manner. At each apex point addi- 
tional diffracted rays are introduced, each of which will travel along the surface of the 
cylinder and will reradiate tangentially. The creeping -wave field for a slotted conducting 
cylinder can be expressed as 
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Figure 6.- Slotted-cylinder creeping-wave fields. 


E c (s, <p) = E^Qj) * -= £ D mh (Q 1 ) D mh (P 1 ) e 0 

V 1 m 

-jk(t 2 +s 2 ) -r* %h(p)<ls 

* E ‘( Q 2) S — ^ Z D -h( Q 2) D nh(P 2 ) e 0 


K only the lowest-order creeping -wave mode is considered, equation (52) reduces to 


E c ( S ,«.E 1 (Q 1 )D oh (Q I )D oh (P 1 )2 


+ E i(Q2) D oh(^2) D oh( p 2) 




and for a circular cylinder simplifies to 


-jk^+r) _ a t 

E c (r>*> = E i( Q l) D oh( Q l) 6 yp e * 

-jk^+r) * 

+ %(Q 2 ) D oh (Q 2 ) e (54) 

since the local radius of curvature is the same at Pp P 2 , Qp and Q 2 . The far-zone 
edge diffracted fields are equated to the creeping -wave fields at the shadow boundary to 
assure continuity of the fields across the boundary between the lit and shadow regions. 
For the geometry of figure 7, the angular variation of the wedge -diffracted fields along 
the shadow boundaries is given by 




e d£' + (! +; ) 

e dHH 


(55) 


(56) 


In other words, the angular variation of the wedge -diffracted fields along the shadow 
boundary is used as the initial value of the angular variations of the creeping waves. 

Thus R* and Rt, serve as the coupling mechanisms between wedge -diffracted and 
creeping waves which assure field continuity along the shadow boundaries. The creeping- 
wave field for a circular cylinder can then be expressed as 


' ik (‘ 1+r ) ‘“oh*! „ /oh 1 ! 


E C < r ,»>° R A yr 6 + yr e 


where 


R a - E i (Q 1 ) D oh(Qi) " E D 


r, + 


(H 


and 


RB=E i( Q 2) D2 h( Q 2) =-^-E D ^-(l4 C ) 


(57) 


(58) 


(59) 
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Figure 7.- Far-zone creeping-wave field coordinates. 

For a noncircular cylinder the diffraction coefficient is a function of the radius of 

of curvature; therefore, the creeping -wave field must be multiplied by D 0 j 1 (P)/d o j 1 (Q) 
to give 


E (s </>) - R ;jN MdS 

E C(=.« AD oh(Ql) Wl 

„ Poh(P 2 ) -i t2ai > h(P)llB 


B D oh( Q 2) \/®2 


( 60 ) 


The wedge approximation of the parallel -plate — cylinder geometry does alter the 
physical boundaries of the antenna structure and the diffraction mechanism of the fields, 
especially in the penumbra region. For the wedge approximation, the fields in the 
( p = ^ directions travel along the boundaries of the wedges. However, in the 
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actual antenna structure some of the energy is trapped on the surface of the cylinder and 
travels along the curved surface reradiating in a tangential direction. Thus the wedge- 
diffracted and creeping -wave fields in the transition region (penumbra) will not satisfy 
the field boundary conditions. However, computations carried out by using such a model 
compare favorably with existing boundary -value solutions and experimental data. There- 
fore, such a model would be a valid approximation of the structure as far as the diffrac- 
tion mechanism of the fields is concerned and equations (55) and (56) would serve as the 
coupling mechanisms between wedge -diffracted and creeping -wave fields. 

CIRCULAR CYLINDER RADIATION 

The wedge diffraction and creeping wave techniques are first applied for the calcu- 
lation of the equatorial radiation patterns of axial and circumferential slots on circular 
conducting cylinders of infinite length. Boundary-value solutions for slots on circular 
cylinders with common feeds exist and are used for comparison. Once the method is 
verified, it is employed for pattern calculations of slots whose field distribution is such 
that modal solutions do not exist and to bodies whose geometric shape does not conform 
to a coordinate system where the wave equation is separable. In addition, the elevation- 
plane pattern for finite-length cylinders is analyzed. It should be pointed out that no 
boundary -value solutions for finite-length cylinders exist, and experimental results are 
used for comparison. The diffraction contributions from the edges of the cylinder and 
the effect of the finite-width aperture to the overall pattern are observed. Pattern calcu- 
lations for more complex geometries such as an elliptical cylinder are carried out in the 
section "Elliptical Cylinder Radiation." 

Equatorial -Plane Pattern of Axial Slots Operating in TEM Mode 

The two-dimensional diffraction geometry for an axial slot mounted on a circular 
cylinder and operating in the TEM mode is shown in figure 8. The junctions formed by 
the walls of the parallel -plate waveguide and the planes tangent to the surface of the 
cylinder at the edge points are represented by a pair of infinite wedges of finite included 
angle, WA = (2 - n)7r. Two regions, lit and shadow, are formed by the imaginary sides 
of the finite wedges as shown in figure 8. 

For 0 g = 90° of figure 3(a), R m (-0 Q ) = 0, R D1R (0 O ) = °> h = W, nj = n g = n, 
then equation (31) for the wedge -diffracted field reduces to 
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LIT REGION 



Figure 8.- Lit and shadow regions for circular cylinder. 

A shift of coordinates from the aperture to the center of the cylinder is convenient, so 
that common coordinates will be used for wedge -diffracted and creeping -wave fields. If 
the far -field approximations 


r Q ~ r - a cos <p cos /3 


(62) 


for phase terms and 


r 


(63) 


for amplitude terms with 




(64) 
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are assumed, the total wedge -diffracted field is expressed as 


E D (r,4>) = 


1 e'j kr 


j k(^sin0+acos0cos|3) - j 

e LV ' 4J |^ 1 (^) + R 2 (0)e-j kwsin f| (65) 


/2^k /r 

The creeping -wave contribution as expressed by equation (57) is given by 


( 66 ) 


The wedge-diffracted and creeping-wave fields, the e"^ kr /\/F factor being sup- 
pressed, are given by 


j k(^- sin 0+acos 0cos /3) - j 


(67) 


and 


■cW) . r a .^ + *) + 


( 68 ) 


In general, the total field is equal to the sum of the two fields 

E t (</>) = E D (<f>) + E c (<f>) (69) 

However, wedge -diffracted and creeping -wave fields do not exist in all regions. To find 
the total field in each region, the appropriate wedge -diffracted and creeping -wave fields 
must be considered as they are tabulated in table I. 

Radiation patterns computed by using the boundary -value solution of Wait (ref. 3) 
and the fields from table I in their respective regions for the diffracted -field solution are 
shown in figures 9 and 10. The second- and higher -order diffractions are approximated 
by wedge diffractions resulting from uniform cylindrical waves. Diffracted waves 
resemble cylindrical waves if they are observed sufficiently far away from any shadow 
boundary that results from their source. The doubly diffracted waves can be adequately 
treated as uniform cylindrical waves for guide widths down to about A/5. However, 
third- and higher-order diffractions are not very accurate when (pg = 90° since they 
are viewed at the shadow boundary formed by their source (ref. 32) and the solution is 
more accurate for wider slots. 
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TABLE I.- WEDGE-DIFFRACTED AND CREEPING-WAVE FIELDS FOR EQUATORIAL-PLANE PATTERN 
OF CIRCULAR CYLINDER IN DIFFERENT REGIONS 


Region I 


Region II 


Region HI 




Region IV 

;+ p < <f> ^ it 


Wedge-diffracted field, 

E D (^) 


jjk(^-5in 0+acos 0CQS/3) 


\J 2nk 

jjk^sin0+ acos <pcospj - 


(R 1 + E 2 e-i kwsin <’) 


\f2irk 

jjk^sin 0+acos (pcosfi^ - 


(R 1 + R 2 e-'i kwsln '*’) 


/2jrk 


-( R i) 


Creeping-wave field, 
E C (<P) 

h 

h 

R A e ’ t l(° l oh +lk ) + R B e' t2 ( < ’ 0h+ik ) 

a (* + F- p ) 


R B e- t ^° feh+lk ) 


-©- 

1 

'XL-' 






R A e-M a oh+i k ) +RBe - t 2(‘WJ k ) 

a (*-M 

a (-»+M 


The solutions which consider only first- and second-order diffractions have a dis- 
continuity in the cf> = ±90° directions because higher-order diffractions are neglected. 
The discontinuity and variation from the boundary -value solution in the 4> = ±90° direc- 
tions are greater for smaller guide widths. The solution which takes into account 
higher-order diffractions eliminates the discontinuity, but it is not very accurate in the 
geometrical shadow boundary region (penumbra) since the assumption of uniform cylin- 
drical waves is not satisfied for third- and higher-order diffractions. 


This analysis is more accurate for large ka cylinders since the decaying expo- 
nents of the surface waves are derived from asymptotic series for large ka. When this 
approach is used, the computations for large ka cylinders do not require any additional 
details from those of small ka. In contrast, the boundary -value -solution computations 
for large ka cylinders are not very convenient since Bessel functions of large order 
and argument are not readily available. Also, the convergence of the radiation pattern 
function becomes poor for large ka, and higher -order terms must be included (about 
2ka terms). From an economic point of view, large ka value computations using the 
boundary -value solution require considerably more computer time compared with the 
corresponding diffracted -field solution (about a factor of 50) since higher -order terms 
must be included for convergence. 


As the value of ka increases, the field in the shadow region 




decreases and the ripples start to appear at larger angles. It is noted that the period of 
the ripples is approximately 180°/ka. A field plot (voltage against angle) for the cylinder 
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with ka = 12 and w/X = 0.2 is shown in figure 11. The variations between the 
boundary -value and wedge -diffracted solutions and the field discontinuity along the 
<p = ±90° directions are more noticeable on a linear scale. The agreement would be 
significantly improved for larger cylinders with wider plate separations as shown for 
the cylinder with ka = 36 and w/X = 0.4 in figure 12. 

Diffracted-field (multiple diffractions) 

Diffracted-field (single-double diffractions) 

Boundary-value 

4 > 

o° 



180 ° 


Figure 11- Field plot of axial infinite slot on circular conducting cylinder (ka = 12, w/X - 0.2) using finite wedges 

(TEM mode). 
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Diffracted-field (multiple diffractions) 
Diffracted-f ield (single-double diffractions) 


Boundary -value 

0 ° 



Figure 12.- Field plot of axial infinite slot on circular conducting cylinder (ka = 36, w/X = 0.4) using finite wedges 

(TEM mode). 

Up to now, the parallel -plate — cylinder geometry has been approximated by a set 
of wedges of finite included angle (n / 2) each formed by a wall of the parallel plate and 
a tangent plane to the cylinder surface at the edge point. For the far field, it would be 
the presence of the conducting wall for the wedge geometry of figure 3 that would force 
the field to be zero at p - it - WA. However, there is no such conducting wall in the far 


field of the geometry of figure 6 and it is only the edge that needs to be considered. 
Sommerfeld (ref. 17) points out that patterns on precise diffraction photographs exhibit 
almost no dependence on the material and shape of the diffraction edge, and a glass sur- 
face with radius of curvature of several meters yields essentially the same diffraction 
fringes as the edge of a razor. Another possible approximation of the parallel -plate — 
cylinder geometry will be a half -plane (wedge with n = 2) instead of a finite wedge 
(n ^ 2). This approximation becomes necessary for the TE^q mode because the wedge- 
diffracted field along the artificial wedge surface separating the lit and shadow regions 
for the finite-wedge approximation is zero. Thus, no creeping waves will be traveling 
around the cylinder surface for the TEjq mode and no fields in the shadow region. Since 
this is not true, the half -plane approximation becomes necessary. 

However, to verify the validity of the approximations - finite wedge (n ^ 2) and 
half -plane (n = 2) - both models are used for the computation of patterns for TEM 
and TE 10 mode slots and compared with boundary -value solutions and experimental 
results. 

It was found, by comparison, that another approximation for the TEM mode slot on 
a circular cylinder was to replace the edge (I) geometry by a half -plane (n = 2) for the 
first-order diffraction and by a finite wedge (n / 2) for second- and higher -order diffrac- 
tions for the 0° to 180° pattern measured in the counterclockwise direction. The edge (2) 
geometry was replaced by a finite wedge for all orders of diffraction. For the 180° 
to 360° pattern, the approximations of edge (I) geometry are valid for edge (2) and vice 
versa. 

The computed results using these approximations are shown in figures 13 and 14 
where they are compared with the boundary -value solution. These approximations 
result in better accuracies for smaller guide widths around the penumbra region than the 
finite-wedge approximation used for figures 9 and 10. However, the finite-wedge approx- 
imation gives better results in the lit region as it should since the fields in the lit region 
are more strongly dependent on the geometry of the edge of the wedge, whereas the far 
fields in the penumbra region would be dependent on the fields on the walls of the wedge. 

It is then concluded that both approximations are valid and a combination of the wedge 
approximations dependent on the location of the observation point would yield optimum 
results. As a first order of approximation, either one of the suggested models may be 
used. A field plot for the cylinder with ka = 12 and w/X = 0.2 cylinder is shown in 
figure 15. The accuracy of this approximation of the penumbra region is better seen in 
this plot when compared with the one in figure 11. 

As was pointed out earlier, many models for the parallel -plate — cylinder geometry 
would give good results as a first-order approximation. Another suggested model for 
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Figure 13.- Radiation patterns of axial infinite slot on smaller circular conducting cylinder using half-plane 

and finite wedges (TEM mode). 




— Diffracted-field (multiple diffractions) 

r Diffracted-field (single-double diffractions) 

Boundary-value 



Figure 15.- Field plot of axial infinite slot on circular conducting cylinder using half-plane and finite wedges (TE/V\ mode). 

the parallel-plate — cylinder geometry would be a set of half -planes (n = 2 for all orders 
of diffractions). Computed results using such a model are shown in figure 16 where rea- 
sonable agreement between the boundary -value and diffracted solutions is indicated. 

Sometimes it may be desirable to have the maximum radiation oriented at a given 
angle other than the zero-degree direction described previously. One way of accom- 
plishing this would be to have the feed of the waveguide mounted in a slanted position. 
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Diffracted-field (multiple diffractions) 
Diffracted-field (single-double diffractions) 
Boundary-value 


— 



180° 


Figure 16.- Radiation pattern of axial infinite slot on circular conducting cylinder using half-planes ITEM mode). 

TEM mode propagation along the parallel plate is assumed but the aperture-field 
distribution will not be uniform in this case. Propagation of this mode in a guide mounted 
on a cylinder in a slanted position has no boundary -value solution. However, the technique 
of diffracted fields can be used to obtain the radiation pattern. The two-dimensional 
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geometry of this mode is shown in figure 17. The modified equation for the diffracted 
field is given by 


E n (0) 


^ j jk acos 0cos/3-|-cos^+e + 0g) "fjr 


Ri(0) 


Dxi " ✓Sk 

4- R 2 (0)e 3khC ° S ^ +e + ^) + Rl (-^) e -j 2kwsin ^ +e ) 


where 




h = 


w 


sin <£ g 


(70) 


(71) 

(72) 



Figure 17.- Slanted feed parallel-plate waveguide geometry. 
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and the regions where the appropriate wedge -diffracted fields are applicable are given 
in the following table: 



The creeping -wave field is the same as that given by equation (68) applicable in the 
regions given in table I. 

The computed radiation patterns for this mode of operation are shown in figure 18 
and a field plot for one cylinder in figure 19. The maximum radiation for the cases 
investigated is oriented at an angle 0 = -e/2. This agrees with the physical reasoning 
by observing the reflection mechanism of the rays from edge ® by the side of wedge 2. 

It can then be concluded that by slanting the feed waveguide at a given angle, it will orient 
the maximum radiation in a direction other than <p = 0°. The discontinuities observed 
in the solution which considers only single-double diffractions occur because third- and 
higher-order diffractions which are neglected become significant. 

Equatorial -Plane Pattern of Circumferential Slots 
Operating in TE^q Mode 

The TE^q mode propagation in a parallel -plate guide can be represented by two 
plane TEM waves reflecting obliquely back and forth between the waveguide walls with 
the electric field parallel to the edge of the walls. The approximation of the parallel- 
plate — cylinder geometry by a set of wedges each formed by a waveguide wall and a 
tangent plane to the cylinder surface at the edge point will not satisfy the tangential 
electric-field boundary conditions. In order to overcome this obstacle, the half -plane 
approximation is used as a model. The boundary -value solution for a circumferential 
slot operating in the TE^q mode and mounted on a circular cylinder exists and is used 
for comparison. 


38 







Diffracted-field (multiple diffractions! 
Diffracted-field (single-double diffractions) 


<p 


0 ° 



270° 


Figure 19.- Field plot of axial infinite slot on circular conducting cylinder with slanted feed (TEM mode). 


The two-dimensional diffraction geometry is shown in figure 20. It was found by 
comparison with the boundary -value solution that one approximation was to replace the 
edge (D geometry by a half -plane (n = 2) for the first-order diffraction and by a finite 
wedge (n / 2) for second- and higher-order diffractions for the 0° to 180° pattern mea- 
sured in the counterclockwise direction. The edge ( 2 ) geometry was replaced by a 
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Observation 

point 



SHADOW REGION 


Figure 20.- Diffraction geometry for a TEio mode circumferential slot. 

finite wedge for all orders of diffraction. For the 180° to 360° pattern, the approxima 
tions of edge (D are valid for edge (2) and vice versa. 

The wedge-diffracted field is given by 


E«(d>) = 


k^sin 0 + acos t^cos^ 


(r, + R,e-i kw sln *) 


(' 



( 74 ) 


R -rW + rM 

K 1 ~ K D1 + K D1 


R 2 = R§2 + R D2 


(75) 


where rJJ}, R^ 1 ], R^, and R^ are given by equations (43), (46), (44), and (47) 
with c£>g = 90°, R D1 (-<£) = 0, and R^ p = 0. The creeping -wave contribution is given 
by 


E(-i(0) — 


-t l ( Q bs + ^ k ) 


+ R p e 


i (°bs + j k ) 


where 
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A “ 



(76) 


(77) 
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as given by equation (73) and 


a. 


os 



(78) 


(79) 


The regions where the appropriate wedge-diffracted and creeping -wave fields are appli- 
cable are identical to those shown in table I for the TEM mode. 

The computed results using this approximation are shown in figure 21 along with 
the boundary -value solution. It is noted that the solution which includes higher -order 
diffractions gives the best results. A noticeable discontinuity is present at the shadow 
boundary when higher-order diffractions are neglected. 

Another approximation of the parallel -plate — cylinder geometry would be a set of 
half -planes (n = 2 for all orders of diffraction) as was true for the TEM mode. Com- 
puted results using this model are shown in figure 22. A good agreement is indicated. 
The discontinuity in the (p = ±90° directions for the single-double diffractions is again 
present since higher -order diffractions are neglected. 
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Diffracted-field (multiple diffractions) 
Diffracted-field (single-double diffractions) 
Boundary -value 



Figure 22 - Radiation pattern of circumferential slot on circular conducting cylinder using half-planes (tEiq mode) 



Elevation -Plane Pattern of Circumferential Slots on 
Finite-Length Cylinders Operating in TEM Mode 


The elevation-plane pattern (F(0 ))^_q for a finite-length conducting cylinder can 
be computed by using wedge -diffraction techniques. The diffracted fields from the ends 
of the cylinder and their contribution to the overall radiation pattern can be readily com- 
puted. Although boundary -value solutions for circular conducting cylinders of infinite 
length exist (refs. 1 and 2), no solution which takes into account finite length is available. 
To check the validity of the technique, experimental results were used for comparison. 

The two-dimensional geometry of a slotted cylinder operating in the TEM mode is 
shown in figure 23. The diffracted field for the TEM mode from edges (D and (D is 
given by 

-i( kr °f J) 





Figure 23.- Diffraction mechanism geometry for elevation-plane pattern. 
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where 


i kw 


i fill sin Or 1 

*>«.)—* °|®i ( s o) + r 2 ( e o) e ”' |kw sto e ° 


( 81 ) 


and and Rg ( ) are given by equations (38) and (34). Diffractions from 

the aperture in the 0 o = ±ir/2 directions will be diffracted by wedges 3 and 4 which in 
turn will be diffracted by wedges 1, 2, 5, and 6. This process continues to higher-orders 
of diffraction. 

The reciprocity principle (see appendix) is applied to calculate the singly diffracted 
field by wedge 3; this leads to 


E D3( r 3 ,e 3) " 


-i( kr 3 + f) 
^/27rkr 3 
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where 


R D3 N ■ b *3- ”3)' 


(83) 


The singly diffracted ray from wedge 3 in the direction of wedge 5 | ' 9 3 = 180°^ will be 
diffracted again and its contribution is 
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D5 ( r 5’ ^5) " 


-j( kr 5 + |) 
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27rkr c 
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(84) 


where 


^ n s)| 


= „(!> 


(85) 


Diffractions from wedge 5 in the directions of wedge 3 = 0°^ and wedge 6 = 270°j 

will cause additional diffractions. If the length of the cylinder is assumed to be several 
wavelengths, the diffractions by wedge 6 are negligible. However, the second-order 
diffractions from wedge 3 are given by 


E D3( r 3> e 3) 




✓ 


R 


27rkr, 


( 2 ) 

D3 




( 86 ) 


46 



I 


where 


(h) =■ R S ( e 5=°)[yB( a * "sj 


(87) 


The same procedure used to calculate the fields diffracted by wedges 3 and 5 can 
be applied to calculate the fields diffracted from wedges 4 and 6 being initiated by the 
aperture diffractions in the 0 O = -tt/2 direction. 

The first-order diffracted fields from wedges 4 and 6 are given by 


where 


E D4( r 4> e 4) = 


-i( kr 4 + |) 
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where 


R D6( e 6) _ R X)4(- 7r ) [yB ( a ’ 27r-0 6 ,n 6)] 


(91) 


and the second -order diffracted field from wedge 4 by 

-j( kr 4 + !) 
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where 


(92) 


R D4 ( 0 4 ) = B ( a > w+0 4* n 4) 


(93) 
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Equations (80) to (93) were all based upon the validity of the reciprocity principle as 
outlined in the appendix. 

Shifting the coordinates from the individual edges to the center of the cylinder and 
assuming the far -field approximations 


r o * r - | cos 9 

rg « r - d cos (y - 9) 

r^ « r + d cos (it - y - 9) ) 

r 5 » r - d cos (it - y - 9) 

rg ~ r + d cos (y - 9) 


(94) 


for phase terms and 


r o sr 3 Kr 4 Kr 5" r 6 sr 


(95) 


for amplitude terms with 


where 


0 o* 0 3 


9 A * 0 , 


9 C ~ 6 



(96) 


(97) 


y= tan -1 - (98) 

cl / 

-j ( kr + j) / 

the diffracted fields, with the e j v/2wkr factor being suppressed, are shown in 

table n. 


It should be pointed out that additional diffractions by each wedge are present and 
can very conveniently be accounted for. However, for cylinders with several wavelength 
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radius and length such diffractions have 
secondary effects and can be neglected. 
It is apparent from figure 23 that the 
two-dimensional geometry used for the 
formulation of the elevation -plane pat- 
tern of a finite -length slotted cylinder 
is identical to the cross-section con- 
figuration of a slotted finite -width 
ground plane. 

The computed results using the 
fields from table II are shown in fig- 
ure 24 where they are compared with 
experimental data. Very good agree- 
ment between theory and experiment 
is indicated. The ripples in the 


TABLE IL- WEDGE-DIFFRACTED FIELDS FOR ELEVATION-PLANE PATTERN 
OF FINITE-LENGTH CYLINDER IN DIFFERENT REGIONS 


Diffraction wedge 

Wedge-diffracted field 

Region 

1 and 2 

jocose 

E D (0)=e 2 [R d (8)] 

-! <e< t 

3 

e dV )(9) = ei kdcos (’'- 0) [R< ) 1 >(e) + Rg>(e)j 

-§<8<„ 

4 

eJ£j 2) (0) = cos(T7-y-g)|li^(0) + R^(8)j 

-»<«<! 

5 

E^(8) = ei kdcos( >'-V-» ) [R^(^ 

o<e<Q 

6 


_ •*£< 0 < o 
2 


- 2. < 9 < §■ region are present because of the diffraction contribution from the edges of 

2i u 

the cylinder, and they become negligible as the length becomes large. The experimental 


results were obtained from slotted ground plane models since it was deduced that their 


two-dimensional geometry is identical to that of a slotted cylinder. 


The beamwidth of the main lobe is decreased as the aperture width is increased. 
The smoothest pattern and the lowest back lobes are obtained for the aperture whose 
width is exactly 1 wavelength since complete cancellation of the fields diffracted from 
wedges 1 and 2 in the 6 = directions occurs. The ripples in the -90° < 9 < 90° 

Cl 

region which are much in evidence in figures 24(a) and (b) do not appear in figures 24(c) 
and (d) because the fields in the 6 - ±~ directions are very weak and any diffractions 

from the edges are negligible. Additional lobes appear as the aperture width is larger 
than 1 wavelength. 


Another formulation for the diffractions from wedges 1 and 2 would be to use con- 
tinuous double diffractions. (See ref. 36.) In this case, doubly diffracted rays from 

wedge 2, R^(0), are neglected in the 0° < 6 < 90° region while doubly diffracted rays 

from wedge 1, R^(0), are neglected in the -90° < 0 < 0° region. Computed results 


using the above formulation are shown in figure 25 where they are compared with experi- 
mental results. Again, very good agreement between theory and experiment is indicated. 
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ELLIPTICAL CYLINDER RADIATION 


Flush-mounted antennas are commonly used on space vehicles, missiles, and air- 
craft. However, their radiation properties usually cannot be predicted analytically 
because the geometric shape of the body to which they are mounted does not conform to 
a coordinate system where the wave equation is separable. The method used to analyze 
the radiation properties of axial slots on circular conducting cylinders can be applied to 
cylindrical -shaped bodies of arbitrary convex cross section. One geometry of wide 
interest is an elliptical cylinder which can be used as an approximate model for the wings 
of an aircraft or the fins of a missile. 

The far-field radiation produced by a slot of arbitrary shape on the surface of an 
elliptical cylinder of infinite length using modal solutions has been carried out by Wait 
(ref. 5). Computations using the modal solution are numerically convenient only for very 
thin-shaped elliptical cylinders, and as the size increases (larger ka and kb) higher - 
order terms must also be included for convergence. However, the analysis using wedge 
diffraction and creeping wave theory can be used for computations of any size and shape 
of an elliptical cylinder. As the physical dimensions increase, the accuracy of the 
obtained data increases because diffraction coefficients Dmh and decay constants 
used in the analysis are derived from asymptotic series of canonical bound ary -value 
problems for X — 0. 

Equatorial -Plane Pattern of Axial Slots Operating in TEM Mode 

The two-dimensional geometry for an axial slot mounted on an elliptical cylinder 
is similar to that for a slot on the circular cylinder and is shown in figure 26 for the 
TEM mode. The junctions formed by the walls of the parallel -plate waveguide and the 
planes tangent to the surface of the cylinder at the edge points are approximated by a 
pair of infinite wedges of finite included angle WA = (2 - n)7r. The wedge -diffracted 
field is given by equation (61) and is repeated here as follows: 



With the use of the far-field approximations 

r Q « r - t q cos <p cos /3 (99) 
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(104) 


X o= ± f 




(105) 


with Rj($) and R 2 (<P) as defined by equations (38) and (34). 

The creeping -wave field is given by equation (60) and is repeated here as follows: 


E °oh(Pl) e ' ik (‘ 1+Sl ) 

c< ’*° A MQT) 


D oh(p 2 ) e -H^) -£ 2 <WW d ‘ l 

BD oh(Q 2 ) f; 6 


+ R 


From the geometry of figure 27, the following far -field approximations may be obtained: 


T ! * t 2 * T 


(106) 


and 


for phase terms and 


5 2 ~ v “ 6 1 


Sj ~ r - r cos 6^ 


s 2 « r + t cos 


(107) 


(108) 


s r s 2 sr 


(109) 


for amplitude terms. Equation (60) then reduces to 
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Figure 27.- Far-zone creeping-wave field coordinates. 




(113) 


R A = ^ E Dp- + (H] * E,(Ql)[l>oh(Qlf 

"b = e d 5 -(| + ?)] = E i(Q 2 )[ D oh(Q 2 f < 114 > 

with Ej-j being obtained from equation (102). 

Expressions for the radius of curvature and differential arc length for an ellipse 
can be obtained by applying the basic geometric definitions of each 


ds = \/(dx) 2 + (dy) 2 

(115) 

[i + (y ’ >5 3/2 

P " |y"l 

(116) 

x = r sin c p 

(117) 

y = t cos cp 

(118) 

Using equations (115), (117), and (118) leads to 


dx = t cos <p d<p + sin <f> dr 

(119) 

dy = -t sin <p d <p + cos <p dr 

(120) 

ds = |r 2 (d(£) 2 + (dr)^] 

(121) 

y » = _( b f sin <t> 

\a/ cos cj> 

(122) 

_ b (a 2 cos 2 (p + b 2 sin 2 cp) ^ 
a 5 cos 2 cp 

(123) 
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Substituting equations (117) and (118) into the equation for an ellipse in rectangular 
coordinates 


leads to 



1 


T = 


ab 


(a 2 cos 2 <p + b 2 sin 2 cp) 


(124) 


(125) 


dT = 


ab 


(a 2 - b 2 ) 


sin 0 cos 0 d ■ 

n /o ' 


(a 2 cos 2 cp + b 2 sin 2 cp)^^ 


(126) 


* 

Substitution of equations (122) and (123) into equation (116) and of equations (125) and (126) 
into equation (121) yields 


P 


l_f a 4 cos 2 cj) + b 4 sin 2 <A ^ 
a k\a 2 cos 2 cp + b 2 sin 2 <p) 


(127) 


ds = ab 


(a 4 cos 2 (j) + b 4 
(a 2 cos 2 0 + b 2 


sin 2 a) 172 
sin 2 ^,) V2 


d(£ 


(128) 


The arc length and complex attenuation factor of the creeping wave on the elliptic surface 
can then be expressed as 


t = r * 12 ds - ab r ^ 12 (a 4 COS 2 <j> + b 4 sin 2 d<p 

^11 ^11 (a 2 cos 2 cp + b 2 sin 2 cp)'* 


(129) 


^ h(p)ds ^-f Jic 


12 
'*11 


[( 


ab d0 


a 4 cos 2 


(p + b 


4 sin 2 


<p){ 


a 2 cos 2 


cp + b 2 sin 2 


<P 


■3 


1/2 


(130) 
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where 


011 


trapping point angle of creeping wave 


012 


exit point angle of creeping wave 


The points of trapping, Qj and Q 2 , are given by equations (104) and (105). To deter- 
mine the points of reradiation, and P 2 , of each creeping wave on the surface of the 

elliptical cylinder, the vector cross product will be utilized (ref. 24). As shown in fig- 
ure 27 for the observation angle 0 with far -field observations being assumed 


h x hj = 0 


(131) 


where 


n = -x cos 0 + y sin 0 = Unit vector perpendicular to line from origin of 

cylinder to observation point 


(132) 


n l = 


A. x •* v 
a 2 b 2 


/ x \ 2 of tangency 

w + U 2 ) 


= Unit vector normal to cylinder surface at point 


(133) 


Applying equation (131) and using the equation for an ellipse in rectangular coordinates 


4 + 4 =i 

a 2 b 2 


lead to, for reradiation point Pp 

x i( p i) = 

y i( p i) = 


-a z cos <fi 


(a 2 cos 2 <p + b 2 sin 2 

b 2 sin 0 

(a 2 cos 2 0 + b 2 sin 2 


- tan-1 


»l( P l) 


(134) 
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and, for reradiation point P 2 , 


i( p 2 ) = 


a 2 cos 


(a 2 cos 2 <p + b 2 sin 2 


H*2) - 


(a 2 cos 2 <p + b 2 sin 2 <p)^^ 


(135) 


The angles /3 and £ are given by 


13 = tan-1 /a w 

\ b \/4a 2 - w 2 


(136) 


£ = tan 


-1 /b w 


^4a 2 - w 2 y 


(137) 


The wedge-diffracted and creeping -wave fields, with the e~j kr //r factor being sup- 
pressed, are given by 


j k(|-sin^ + T 0 cos^cos/3') 

E ° W = 7fcr L J ( R ! + R 2 e- ikWS ‘ n ") 


(138) 


E n (0) = P A e e 


-jktj -jkt 2 -£ 2o bhM ds 


+ P^e " e 


(139) 


tl = r^ 12 ds 


(140) 


■£ 


'022 

t 2 = \ ds 


0 


21 


In general, the total field is given by 


(141) 


E t (0) = E D (<f>) + E c (0) 


(142) 
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where the appropriate fields must be considered in each region as shown in table m. 


TABLE HI.- WEDGE-DIFFRACTED AND CREEPING-WAVE FIELDS FOR EQUATORIAL -PLANE PATTERN 
OF ELLIPTICAL CYLINDER IN DIFFERENT REGIONS 


Region I 

0<<M§- C 


Region H 


Region m 


Region IV 

If+CcMir 


Wedge-diffracted field, 

E d (0) 


j|k^|-sin0+T o cos^cos^ - | 


jjk^sin 0+t o cos 0cos/3^ - |j 


jjk^sin 0+ r 0 cos 0cos ^ 




(R 1 + R 2 e-^ sin ^ 
(r x + R 2 e“^ kw sin *) 
>l) 


Creeping -wave field, 
E c (0) 


-JW 1 -J* lo bhW ds -jkt 2 -j’‘ 2 °bhW ds 

P A e L e u +P B e e 

-jkt 2e -l‘ 2 “ohW ds 

-iktj \f“oi>> ds 
e 

-jktj -j’‘ 1 “ohW ds -jkt 2 -Io 2 °ohW ds 


*11 

*21 

*12 

*22 

*n = <* 

-©• 

CO 

11 

1 

0 12 = 2ff + 

*22 = -0 


*21= -(»‘*l) 


*22 = 


■©- 

CO 

11 

=4 


*22 = _/3 

*n = d 

*21= - *l) 

*12 - *1 

*22 ' - s 


J 


Radiation patterns computed by using the fields from table m in their respective 
regions are shown in figures 28 to 30. Experimental models were constructed and the 
obtained data were used for comparison since boundary -value solutions are numerically 
convenient only for very thin elliptical cylinder. Harrington (ref. 4) points out that the 
equatorial radiation pattern (0 = 90° plane) for a circular cylinder with finite axial slot 
is identical to that for a circular cylinder with an infinite axial slot. It is assumed that 
such a relation holds for an elliptical cylinder, and the experimental results were obtained 
by using finite axial slots. The experimental models were also of finite length of about 
30 wavelengths; diffractions from the edges of the cylinders should then be negligible. 
Dimension construction accuracies were within ±X/25 and surface irregularities within 
±A/50. All measurements were performed at a frequency of about 10 GHz. 

The agreement between computed and experimental patterns is good particularly in 
the forward hemisphere. Errors of the order of 2 to 3 dB are obtained in the penumbra 
region and they indicate that further improvement could be made on the approximations of 
the parallel -plate — cylinder geometry by a wedge and on those implicit in equations (113) 
and (114). The use of the half -plane wall for certain diffraction processes would yield 
improved results near the penumbra region as has been discussed previously for the cir- 
cular cylinder. In general, the errors in the shadow region are of the same order of 
magnitude, and they are quite acceptable for determining possible interference and noise 
that might be introduced from other systems through these back lobes. 
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Figure 28.- Radiation patterns of axial infinite slot on elliptical conducting cylinder using finite wedges (TEM mode) 






I 


For elliptical cylinders with constant major axis but decreasing minor axis (thinner 
shape ellipse), the rate of attenuation of the creeping waves, whose main contribution is 
in the shadow region, is larger. However, the lit-region radiation, whose primary source 
is the wedge -diffracted fields, is essentially unaffected as shown in figures 28(a) to (c). 
The creeping -wave fields are strongly influenced by the radius of curvature which is 
rapidly changing whereas the wedge -diffracted fields depend on the included angle of the 
wedges WA which remains essentially constant. 

For the small elliptical cylinder (ka = 10.666, kb = 8) of figure 29(a), slightly 
larger variations and asymmetries between experiment and theory are noted. For 
±1.0 percent changes in frequency (figs. 29(b) and (c)), the variations and asymmetries 
are even greater and it is believed that the physical construction tolerances of the ellip- 
tical cross section are critical for the measurements on this size body. For a thin 
elliptical cylinder the radius of curvature of the surface near the major axis is small 
and any diffractions from that region would not be very accurate since the asymptotic 
series assumption of large physical dimensions is not well satisfied. 

In figure 30, the experimental pattern of a thin elliptical cylinder (ka = 40, kb = 8) 
is compared with a theoretical curve of a finite -size ground plane obtained by using wedge 
diffraction techniques. A good agreement is noted in the 90° < 0 < -90° region. How- 
ever, in the 90° > <p > -90° region the ripples predicted by the computed curve of the 
ground plane do not appear in the experimental curve for the thin elliptical cylinder. 

The surface of the cylinder is smooth and has no sharp edges to contribute diffractions 
in the lit region for the ripples to appear. The fact that edge-diffracted fields do disturb 
the pattern is interesting to the antenna designer who would certainly prefer the smooth 
pattern. Cylindrical caps on the sides of the ground plane would reduce the ripples and 
smooth the pattern. 

To check the validity of the half-plane (n = 2) approximation for the parallel -plate — 
cylinder geometry instead of a finite wedge (n ^ 2), equations (138) and (139) were used 
with modified n's and the computed results were compared with experimental results. 

It was found by comparison, like for the circular cylinder, that the best approximation 
was to replace the edge (1) geometry by a half -plane (n = 2) for the first-order diffraction 
and by a finite wedge (n ^ 2) for second- and higher-order diffractions for the 0° to 180° 
pattern measured in the counterclockwise direction. The edge (2) geometry was replaced 
by a finite wedge for all orders of diffraction. For the 180° to 360° pattern, the approxi- 
mations of edge CD geometry are valid for edge (2) and vice versa. 
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The computed results using the half -plane and finite-wedge approximation are 
shown in figure 31, where they are compared with the experimental curves. It is 
seen that these approximations result in better accuracies in the penumbra region 
than the finite wedge approximation but larger variations exist in all other regions. 

As concluded previously, the solution used then depends on the region of space of 
interest. It should be noted that computed results for single-double diffractions and 
multiple diffractions shown in figures 28, 29, and 31 are almost the same, and it is 
rather difficult to distinguish any variations between the two. 

Equatorial -Plane Pattern of Circumferential Slots 
Operating in TEjq Mode 

As was explained previously, the TEjq mode is represented by two plane TEM 
mode waves reflecting obliquely back and forth between the waveguide walls with the 
electric field being parallel to the edge of the walls. The finite -wedge approximation 
for the parallel -plate — cylinder geometry will not satisfy the boundary conditions and 
the half-plane approximation will be used. 

For a circumferential slot operating in the TE-^q mode and mounted on a circular 
cylinder it was found that the best approximation was to replace edge © by a half- 
plane (n = 2) for the first-order diffraction and by a finite wedge (n ^ 2) for the second- 
and higher-order diffractions for the 0° to 180° pattern measured in the counterclock- 
wise direction. The edge © geometry was replaced by a finite wedge for all orders 
of diffraction. For the 180° to 360° pattern, the approximations of edge (1) geometry 
are valid for edge © and vice versa. 

This approximation is also used for the elliptical cylinder. The wedge -diffracted 
field is given by equations (73) to (75) with n modified as explained previously. The 
creeping -wave field is given by 


ch 

"3^1 Jn 
E c (0) = P A e 1 e U 


% s (p)ds 


+ P B e 


— jktr, 


_f*2 

. 4 > 


<* os (p)ds 


(143) 
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where 





9 4 9 

cos <p + t>* sin^ 


ab d0 

4>){a? cos^ (f> + sin^ 



(144) 


The regions where the appropriate wedge -diffracted and creeping -wave fields are appli- 
cable are identical to those shown in table HI for the TEM mode. 

The computed results for the TE^q mode slot on an elliptical cylinder are shown in 
figure 32. As for the circular cylinder, the single-double diffraction solution has a dis- 
continuity in the <p = ±90° directions because higher -order diffractions are neglected. 
The thin -walled guide approximation (n = 2 for all orders of diffraction) could also be 
used as was demonstrated for the circular cylinder. Good results would also be obtained 
in this case. 


CONCLUDING REMARKS 

Wedge diffraction and creeping wave theories have been used to analyze the radia- 
tion properties of slots on cylindrical bodies. By using this technique, the computed 
patterns for circular and elliptical cylinders have been favorably compared with existing 
modal solutions and experimental results. 

The concept of a hybrid solution utilizing wedge diffraction and creeping wave 
theories has been established in solving aperture antenna problems on cylindrical bodies 
that have not been considered previously from this point of view. The first attempt was 
to apply this technique to a circular cylinder with a known boundary -value solution for 
comparison of the results. It was then extended to a more complicated geometry such 
as the elliptical cylinder where the computations using the boundary -value solution are 
not very convenient but where experimental data are available. Radiation properties of 
finite-size cylinders were also analyzed and the contribution to the overall pattern from 
the edges was taken into account. 

Once a good understanding of the radiation mechanism is established, the technique 
can now be extended to other more complex geometries to include slots on spheres, pro- 
late spherioids, and ogives. Problems involving dipoles on the surface of such geo- 
metric shapes should also be considered. Experimental data and modern computational 
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methods such as the point- matching boundary -value technique should be used in combi- 
nation with the geometrical theory of diffraction to obtain improved values and gain 
better understanding of the radiation mechanisms involved. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 13, 1969. 
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APPENDIX 


DIFFRACTION BY A PERFECTLY CONDUCTING WEDGE 


Consider a two-dimensional wedge and a line source to be situated in space as 
shown in figure 33. The faces of the wedge are formed by two semi-infinite intersecting 
planes. The infinitely long line source is parallel to the edge of the wedge, and its posi- 
tion is described by the coordinate (r 0 ,i[/ o y The typical field point is denoted by 
The line source is assumed to have unit strength and time dependence of the form e^ . 
The Green's function for this radiating system can be written as (ref. 19) 


OO 



m=0 


e m 'Jm(kr) (kr Q ) cos n (4 / ~ ^o) 


n 


± cos ^ + p o ) 


(145) 


if r Q > r. For the case r Q < r, r and r 0 are interchanged. In this expression e m 
is the Neumann number which is equal to 1 if m is 0; otherwise it is equal to 2. The 
plus sign between the two cosine terms is used if the boundary condition is of the homo- 


geneous Neumann type 



on both faces of the wedge). 


For the homogeneous 


Dirichlet boundary condition (G = 0 on both faces of the wedge), the minus sign is used. 


Field point 



Figure 33.- Line source and diffracting wedge. 
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This convergent series is an exact solution to the time harmonic, inhomogeneous wave 
equation for the problem of a radiating line source and wedge embedded in a linear, iso- 
tropic, homogeneous, lossless medium. 

Many times it is necessary to determine the total radiation field when the line 
source is far removed from the vertex of the wedge. In such cases, equation (145) can be 
simplified by replacing the Hankel function by the first term of its asymptotic expansion, 
that is, by the relation 



-j 

e 



1 m 7[\ 
°"4" n 2 ) 


(146) 


This substitution reduces G to the form 



e m Jm( kr ) e 
n 


+ 1 ™ 1 
J n 2 


COS ^(l// - ip o ) 


± COS ^(l// + ^o) 


(147) 


where 


i v i— 1 r 

E = E d + E g = - 2, e m JmW e H 2 cos " ^ Q ) ± cos + i(/ Q ) 
™_n n L 


m=0 


(148) 


with Ep> and Eq being the total diffracted and geometrical optics fields, respectively, 
and E, the series form of the Green's function describing the total field created by the 
diffraction of a plane wave by a wedge. 

An asymptotic expansion for E in inverse powers of kr is very useful for com- 
putational purposes, because of the slow convergence of equation (148) for large values 
of kr. In order to derive an asymptotic expression for E by the standard method of 
steepest descents, it must first be transformed into an integral or integrals of the form 


F(z)e kr dz 


(149) 
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and then evaluated for large kr by means of the method of steepest descents. Such a 
procedure leads to the expression for the diffracted field for a plane wave incident on a 
wedge of included angle (2 - n)n shown in figure 34 to be (ref. 19) 


ED(r,^ 0 > n ) = v B( r > V'-V'o, n ) ± V B (r, ^+^ 0 , n) (150) 

where the plus sign applies for the polarization of the electric field perpendicular to the 
edge 



= 0 


and the minus sign applies for polarization parallel to the edge 


^wedge - 0 



Figure 34.- Diffraction by a wedge of included angle (2 - n )ir. 
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The wedge -diffraction function V B (r, \f/+\f/ 0 , n) for a plane wave has also been deter- 
mined by Pauli in reference 18. 

The form of the diffraction function is defined as 


V B (r, <W 0 , n) = I.^(r, xfy+xp Q , n) + I +7r (r, \p+xf/ 0 , n) 


( 151 ) 


where 


W r > M 0 > n > 


^ e 


-j(kr+|) 


]n\j2ir 


\fg cot 


+ * Q ) jkrg n 


2n 


\ e -* T dr + Higher -order terms 

•J\fkrg 

(152) 


g = 1 + 


cos (xj/ ± ^oj - 2n7rN 


(153) 


and N is a positive or negative integer or zero which most nearly satisfies the equation 


2 n 7 rN - (x// ± = -n 

2nwN - (if/ ± iI/qJ = +7r 


(for I.*) 
(for I +7r ) 


(154) 


Equation (152) contains the leading terms plus higher -order terms which are negli- 
gible for large values of kr. For large values of krg, equation (151) reduces to the 
form presented by Pauli (ref. 18) and given by 


V B (r, xf/+\p 0 , n) 


4 »?) 


5-“S 


COS I - cos 


(155) 


The diffracted field of equation (155) is that from which the asymptotic diffraction coeffi- 
cients of the geometrical theory of diffraction are obtained. (See ref. 10.) This expres- 
sion is not valid in the shadow boundary because g = 0 there. The geometrical optics 
field Eq is given by equations (5) to (7). 

The solution for cylindrical wave diffraction at large distances from the edge can 
be determined by the use of the principle of reciprocity together with the solution for 
plane wave diffraction. (See ref. 6.) For plane wave incidence the diffracted field at 
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observation point P of cylindrical coordinates (r,i //) as shown in figure 35(a) is given 
by Vg(r, n) of equation (151). Now consider the situation in figure 35(b) in which 

the wedge is illuminated by a cylindrical wave with its source at (x Q , By reciprocity 
the diffracted field in the direction | is equal to the diffracted field Vg which is 

located at the point (r=x 0 , with a plane wave incident from the direction ip Q = 
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